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We investigate tlie dynamics of Eulerian walkers as a model of self-organized criticality. The 
evolution of the system is subdivided into characteristic periods which can be seen as avalanches. The 
structure of avalanches is described and the critical exponent in the distribution of first avalanches 
T = 2 is determined. We also study a mean square displacement of Eulerian walkers and obtain a 
simple diffusion law in the critical state. The evolution of underlying medium from a random state 
to the critical one is also described. 



I. INTRODUCTION. 

To illustrate the phenomenon of the self-organized crit- 
icality (SOC) § a wide range of cellular automata such 
as sand piles, rice piles and forest fires, have been pro- 
posed . They assume a system consisting of a large 
amount of elements. The energy being randomly added 
to the system is redistributed among the degrees of free- 
dom by a kind of nonlinear diffusion. This is realized by 
means of avalanche-like processes which carry the added 
energy out of the system. As a rule, the system sponta- 
neously evolves towards the critical state free of charac- 
teristic length and time scale. In this state, probabilistic 
distributions of quantities, characterizing the statistical 
ensemble, exhibit the power law behavior. 

What features of the SOC dynamics are responsible 
for the existance of a dynamical attractor in complex 
systems? What are the origins of the scaling and self- 
similarity in the stationary state that appeares? To an- 
swer these questions, one evidently needs investigating 
nonlinear diffusion in the SOC models and studying the 
structure of avalanches. Due to a rather complicated 
dynamics of most SOC models, the description of their 
evolution is a very difhcult problem. Up to now, the 
most analytically tractable model has been the Abelian 
sandpile model (ASM) Due to its simple algebraic 
structure, the detailed description of the SOC state of 
ASM has been given, and some critical exponents have 
been found |-§. 

Recently, a new model has been proposed which is 
called the Eulerian walkers model (EWM) It also 
possesses Abelian properties and exhibits SOC. While 
EWM has many common features with ASM, its dy- 
namics is much more transparent. This model is formu- 
lated as a deterministic walk interacting with a change- 
able medium. The motion of a particle is affected by 
the medium, and in its turn affects the medium inducing 
strong correlations in the system. If the walk occurs in 
a closed system, it continues infinitely long and eventu- 
ally gets self-organized into Eulerian trails. If a system is 
open, the particles can leave the system and new parti- 
cles drop time after time. In this case, the system evolves 
to the critical state similar to that in ASM. By analogy 
with ASM, the avalanches in EWM have been introduced 



[|10| as periods of reconstruction of recurrent states, after 
they have broken by an added particle. 

In this article, we use the properties of Eulerian trails 
to describe the evolution of the system in detail. In this 
way, we find the critical exponent in the distribution of 
first avalanches. Using the Green function of EWM in- 
troduced in 1^], we obtain the critical exponent charac- 
terizing spread of particles with time. 



II. ALGEBRAIC PROPERTIES OF EULERIAN 
WALKERS MODEL. 

The Eulerian walkers model is defined as follows. Con- 
sider an arbitrary connected graph G consisting of N 
sites. Each site of G is associated with an arrow which 
is directed along one of the incident bonds. The arrow 
directions at the site i are specified by the integers rii, 
(1 < Tii < Ti) where is the number of nearest neighbors 
of the site i. The set {ui} gives a complete description 
of the medium. Starting with an arbitrary arrow config- 
uration one drops the particle to a site of G chosen at 
random. At each time step: 

(i) the particle arriving at a site i changes the arrow di- 
rection from rii to ni + l{modTi) 

(a) the particle moves one step along the new arrow di- 
rection from i to the neighboring site i' . 

Having no endpoints on G, the particle continues to 
walk infinitely long. Due to a finite number of possi- 
ble states of the system, it eventually settles into the 
Poincare cycle. For most dynamic systems the recurrence 
time of this cycle grows exponentially with N. It has been 
shown in ||^ that for the Eulerian walk the Poincare cycle 
is squeezed to the Eulerian trail with the recurrence 
time of an order of N. During the Eulerian trail the par- 
ticle passes all bonds of the graph exactly once in each 
direction. 

Let G be an open graph. It means that one auxiliary 
site is introduced which is called a sink. The subset of 
sites of G connected with the sink forms an open bound- 
ary. The sink does not have an arrow and the particle 
reaching the sink leaves the system. After that, a new 
particle is dropped to the a site of G chosen at random. 
Since on the closed graph the particle visits all sites dur- 
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ing the walk, at the open one it ever reaches the sink. A 
set {C} of configurations C = {rij} which remains on G 
after the particle left G for the sink, is the set of stable 
configurations. The operator can be introduced 

a^C = C (1) 

which describes the resulting transformation caused by 
dropping the particle to the site i. As usual in the the- 
ory of Markov chains, the set {C} may be divided into 
two subsets. The first subset denoted by {R} includes 
those configurations which can be obtained from an ar- 
bitrary configuration by a sequential action by the oper- 
ators Qi. It follows from the definition that the subset 
{R} is closed under a multiple action by the operators 
Qi. Once the system gets into {R}, it never gets out 
under subsequent evolution. All nonrecurrent configura- 
tions are called transient and form the subset {T} which 
is the complement to the set {R}- By definition any re- 
current configuration C € {R} may be reached from any 
another C G {R} by a subsequent action of the opera- 
tors a;. Since this is valid for C — C too, the identity 
operator acting in {R} exists. In addition, the operators 
tti have the following properties: 

1. For arbitrary sites i and j and for any configuration 
of arrows C 

QiGjC = ajQiC (2) 

2. For any recurrent configuration C £ {R}, there 
exists a unique 

(a-'C) € {R} 

such that 

a,(a-i)C = C (3) 

The proof of these statements is similar to ones for the 
avalanche operators in ASM and is given in §. Thus, 
the operators Oj acting in the set of recurrent configu- 
rations {R} form the Abelian group. The addition of n 
particles to site i gives the same effect as the addition of 
one particle to each of neighbors of i. It returns the 
arrow outgoing from i to the former position and initi- 
ates the motion of one particle to each neighboring site. 
In the operator form this looks like 

fc=i 

where jk are the neighbors of the site i. Introducing the 
discrete Laplacian on G as 

{n , i = j 

— and j are connected by bond (5) 
, otherwise. 

and using Eq.(^), one can write the identity operator as 



Since all recurrent configurations can be obtained from 
an arbitrary one by a successive action by operators a^, 
one can represent any C G {R} in the form 

C^l[{a,rC*. (7) 

ieG 

The A^-dimensional vector n labels all possible recur- 
rent configurations. Eq. shows that two vectors n 
and n' label the same configuration if the difference be- 
tween them is mjAij where mj are integers. The TV- 
dimensional space {n} has a periodic structure with an 
elementary cell of the form of a hyper-parallepiped with 
base edges = (Aa, Ai2, AiA?). Thus, the number of 
non-equivalent recurrent configurations is 

N = detA, (8) 

which is Kirchhoff's formula |l^] for spanning trees and 
Dhar's formula for ASM [Q. The correspondence to ASM 
is not surprising. The algebra of the operators com- 
pletely coincides with that of avalanche operators of the 
Abelian sandpile model [Q. Moreover, the identity oper- 
ator has the same form for both the models . This is 
the reason why the numbers of recurrent configurations 
coincide. 

Continuing the analogy between EWM and sandpiles, 
one can find the expected number Gij of full rotations of 
the arrow at site j, due to the particle dropped at i [Q. 
During the walk, the expected number of steps from j is 
AjjGij whereas —J2k^j GikAkj is the average fiux into 
j. Equating both the fluxes, one gets 

^G.feA,, (9) 

or 

G., = [A-%. (10) 

The expected number of full rotations of the arrow is 
equal to the number of entries into the site j divided by 
Tj. On the other hand, the number of visits of a site for 
the random walk is also Green function of the Laplace 
equation. Thus, we have a surprising fact that the num- 
ber of visits of the site by the particle for the determinis- 
tic motion in EWM coincides with that for the ordinary 
random walk. 

The direct correspondence between spanning trees, re- 
current configurations of EWM and Eulerian trails can 
be established in the following way. The walking par- 
ticle leaves each site along an arrow after turning the 
arrow. Therefore, the trajectory of the particle is a trace 
of arrows. If G is an open graph, all trajectories end at 
the sink and never form loops. The corresponding arrow 
configurations are acyclic ones. 
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Given an acyclic arrow configuration, we can construct 
a unique spanning tree rooted in the sink and vice versa. 
Indeed, all bonds along which the arrows are directed 
form the spanning tree. Conversely, if we have a span- 
ning tree rooted in the sink we can obtain the acyclic 
arrow configuration by pointing the arrow from each site 
along the path leading to the sink. This correspondence 
allows us to identify the acyclic arrow configurations and 
spanning trees. Further, we do not distinguish between 
them. Saying spanning tree we mean both the spanning 
tree and its arrow representation. 

If G is the closed graph, the particle settles into the 
Eulerian trail during which it passes each bond exactly 
once in each direction. Let the particle, which has al- 
ready visited all sites, arrive at a site i at some moment. 
If we now remove the arrow from z, we obtain the acyclic 
arrow configuration where arrow paths from any site end 
in i. This defines the spanning tree rooted in the site 
of the current particle location. Therefore, given a site 
i, each Eulerian trail on a closed graph has the unique 
spanning tree representation. 



III. AVALANCHE DYNAMICS. 

The particle added to the recurrent configuration of 
ASM may induce successive topplings of sites called the 
avalanche. At the initial moment, it destroys the recur- 
rent configuration and the system leaves the critical state. 
After the avalanche stops, the recurrent configuration is 
restored again. Thus, the avalanche in ASM may be de- 
fined as a period of reconstruction of the recurrent state. 
This definition may be directly applied to EWM. 

We start with a recurrent state of EWM. The cor- 
responding arrow configuration forms a spanning tree. 
Once a particle is dropped, it may destroy the spanning 
tree by closing a loop of arrows. During the evolution, 
one loop can be transformed into another. When all loops 
disappear, the spanning tree is restored. The interval of 
existance of the loop can be called the avalanche of cyclic- 
ity or simply avalanche. The loops may be created and 
destroyed several times during the motion of one parti- 
cle. Therefore, unlike ASM, an addition of one particle 
may initiate several avalanches in the system. When a 
particle comes to the sink, it always directs the arrow 
to the sink thus restoring the spanning tree. Therefore, 
when the particle leaves the system the avalanche always 
ends and recurrent state is restored. All motions of the 
particle represent successive transitions from one recur- 
rent state to another through avalanches. To study the 
evolution of the system, the structure of the avalanche 
should be considered in detail. 

Consider the Eulerian walk on the square lattice C of 
size L X L with open boundary conditions. Each bound- 
ary site is connected to the sink by one bond on the edge 
and by two bonds at the corners of £. The rule of arrow 
rotations is the same for all sites. If we denote the bonds 



outgoing from a site i by A'', E, S, W, the rule of rotations 
is TV — > i? — > iS — > W. In other words, when the particle 
arrives at a site, the arrow outgoing from this site turns 
to the next bond clockwise. Due to a topological rea- 
son, this rule leads to a simple structure of avalanches, 
namely to compactness of clusters of sites visited by the 
particle. 

Let the particle be dropped to a recurrent configura- 
tion which is a spanning tree. At some step the first 
loop is created. The arrows can form loops of two kinds: 
clockwise and anti-clockwise. The loop is clockwise if 
tracing the loop along arrows leaves the interior of the 
loop on the right and anti-clockwise otherwise. It is easy 
to see that due to the clockwise rule of rotations, only 
clockwise loops may be created from recurrent states. 
Indeed, the anti-clockwise loop arises when the arrow, 
which closes this loop, is directed at the previous time 
step into the area bounded by the loop. The arrow path 
beginning from this arrow could not leave the area of 
the loop without intersections with the loop. This means 
that before this loop was closed, another loop existed, 
which contradicts the assumption that we start with a 
spanning tree. 

Consider the evolution after the closing of a clockwise 
loop at the spanning tree. Denote by ij the arrow if it is 
pointed from site i to site j. Analogously, we denote by 
1112*3... the arrow path if the arrow from site ii is pointed 
to site ^2, the arrow from i2 is pointed to 13 and so on. 
Let a spanning tree exist at the time step {t — 1), while 
at the step t, the particle arrived at the site ii changes 
the arrow direction from iiIq to iii2 and the clockwise 
loop — 111213. ..Zn*i appears (Figj^a). Now, we can 
prove the following: 

Proposition 1: The particle does not leave the area of 
the loop 0+ and the spanning tree cannot be restored un- 
til all arrows inside the loop area make the full rotation 
and the arrows belonging the loop itself change the di- 
rection to anti-clockwise forming the anti-clockwise loop 
= iii„...i2ii. At the last step when 0~ appears, the 
particle arrives at 12 and at the next step the arrow at i2 
rotates out of the loop area and the spanning tree may 
be restored (Fig. |l|b). 

Proof: Consider EWM on the auxiliary graph Q, 
which is a part of the square lattice bounded by the 
loop with closed boundary conditions. The closed 
boundary means that all bonds that link boundary sites 
ii, Z2, ia, with the sites of the lattice outside the loop 
area are removed. The rules of rotations are modified in 
such a way that an arrow skips deleted bonds. We con- 
sider the Eulerian trail at Q starting from the site Z2. 
At the initial moment, the arrow configuration at Q dif- 
fers from that on the lattice C only by orientation of the 
loop: instead of the clockwise loop = 111213. ..i„ii on 
C, we have the anti-clockwise loop 0~ = iii„...i2ii on 
G (fig.|l|c). Starting from the first step, (n — 1) succes- 
sive steps reverse 0~ into 0+ and the particle arrives 
at ii (fig.|l|d). Notice that the initial arrow configura- 
tion on Q corresponds to that described in the previous 
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section, when the particle has already settled into the 
Eulerian trail on the closed graph. Indeed, at the first 
moment, all arrows except the arrow at the current par- 
ticle location site form the spanning tree rooted in this 
site. Hence, the subsequent evolution leads again to the 
loop via full rotation of arrows at all internal sites 
(fig.|l|e). On the other hand, this part of evolution of the 
graph Q coincides with one on the original lattice C since 
the moment when the loop 0+ is closed (Fig. |l|a) up to 
the moment when it is changed by 0~ (Fig. ||b). At the 
last step 12*1 rotates out of the loop area and the loop 
may be broken. Before this moment the loop exists per- 
manently as during Eulerian trail one loop always exists. 
The proposition is proved. 

Generally, the avalanche docs not necessarily end af- 
ter that. Two situations are possible. At the last step, 
the arrow at ?2 turns outside the anti-clockwise loop 
*2*i ~* ^2^2- If is connected to the sink through the ar- 
row path, the spanning tree is restored and the avalanche 
is finished. This is the case of a one-loop avalanche. In 
the other case, the arrow path from i'2 goes to 12, i.e. i'^ is 
the predecessor of with respect to the sink. Then, one 
more loop is closed and the avalanche continues. This is 
a two-loop avalanche. The second loop relaxes like the 
first one. When the second loop is reversed, the span- 
ning tree is always restored because at the last step the 
particle arrives at zq which was connected to the sink by 
an arrow path before the avalanche started. 

Several consequences may be obtained from the picture 
described. During the avalanche the particle visits sites 
inside the loop four times, sites of the edge two times, 
one time at the corner ^ and three times at the corner 
Then, the duration of relaxation of a loop is given 
by the formula 



T = (4s + 2p - 4) + 1 



(11) 



where s is the number of inner sites, and p is the perime- 
ter of the loop. As the avalanches may consist of one or 
two loops, the duration of avalanches may be equal to 



Ti = 2/c + 1 



(12) 



where fc = 0, 1, 2, ... This explains the double distribution 
of durations of avalanches (Fig. ||) obtained in ||l^ . Also 
we can find the critical exponent of the duration distribu- 
tion for the first avalanche. In the thermodynamic limit, 
the duration of avalanches grows as the area of the loop. 
It has been shown in that the probability to get a 
loop of the size s when a bond is added to the spanning 
tree at random is equal to 



(13) 



The distinction of our case from this one is that the loop 
is closed by turning the single arrow that was connected 
to the sink through an arrow path before the turn. Hence, 



the distribution (13) should be divided by the perime- 
ter of the loop. Taking into account that that perime- 
ter scales with the linear size as a fractal dimension of 
a chemical path p ^ and that the loop is compact 
s ~ r^, we obtain 



(14) 



Thus, for the first avalanches the critical exponent of the 
distribution of duration is r = 2. The one- and two-loop 
avalanches differ only in a local structure of the spanning 
tree at the site of closing the loop. Therefore, the critical 
exponents are the same for both the distributions. This 
result is in excellent agreement with numerical simula- 
tions presented in Fig. where we have considered the 
EWM on the square lattice of linear size L — 400 with 
open boundary conditions. 

The evaluation of the exponent of the first avalanche 
size distribution in EWM is similar to one for the distri- 
bution of sizes of erased loops in the loop-erased walks, 
which was studied in [ p^ . The same exponent r = 2 was 
obtained. 

The result t = 2 is valid only for the first avalanches for 
their independence of each other. The analytical deriva- 
tion of T for arbitrary avalanches is a more difficult prob- 
lem due to correlations between subsequent avalanches 
appearing during the evolution of one particle. 



IV. PROPAGATION OF EULERIAN WALKERS. 

Besides the evolution of the system as a whole, we can 
describe the motion of the particle itself. Consider the 
particle dropped on the lattice with a spanning tree. We 
call the site i a predecessor of j if the arrow path comes 
from i to J. Since the particle motion is traced by a path 
of arrows, all visited sites are predecessors of the site of 
a current particle location. If the particle arrives at the 
site which is its predecessor, the loop is closed. Thus, the 
particle can visit the sites that have already been visited 
only during an avalanche. 

We subdivide the motion of the particle into the fol- 
lowing stages. The first stage coincides with the first 
avalanche. At the moment it finishes, the avalanche area 
remains bounded by the anti-clockwise loop opened at 
the bond connecting two sites where it begins and where 
it ends. Further, moving on the lattice, the particle can- 
not enter the area of the first avalanche. New avalanches 
appear beyond the first one being attached to its bound- 
ary and tending to go clockwise around it. Eventually, 
the particle creates a loop enclosing the area of the first 
avalanche. When the avalanche corresponding to this 
loop ends, the second stage of the evolution finishes. At 
this moment, we have the cluster of visited sites which 
consists of the area of the first avalanche, where each in- 
ner site is visited eight times, surrounded by the clusters 
of subsequent avalanches, where all sites are visited four 
times (Fig. |). 
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The further behavior of the system is similar. If at 
some evolution stage we have a cluster of visited sites, at 
the next stage all sites of this cluster will be visited four 
more times and some new area will be added to the clus- 
ter of visited sites. After each evolution stage finishes, 
the cluster of visited sites is compact because it consists 
of compactly situated avalanche clusters. 

Thus, we obtain the system of compact clusters where 
the sites are visited 4iV, = 1, 2, ... times. The clusters 
are strictly embedded one into another with a growing 
number of visits like Grassberger-Manna clusters in ASM 

Using this picture, we can find time dependence of the 
mean square displacement of the particle in the critical 
state. The number of visits {N{R)) of a site separated 
from the origin by the distance R is given by the Green 
hmction of the Laplace equation Eq.(^). When |r — r'| 
tends to the lattice size, G(r, r') decays as log(L/|r — r'|), 
so we can write 



dN{R) 1 
dR R' 



(15) 



On the other hand, the time T required for the parti- 
cle to visit four times all the sites of the compact cluster, 
is of an order of its size R^. Then, the velocity of the 
growth is 



dN 



1 

i?2 



(16) 



Using (^ and (|l^) and the property of the embedded 
clusters, we obtain the mean square displacement 



< R-" >r 



(17) 



that is the diffusion law of the simple random walk. 

In the transient state, we have no the spanning tree 
representing the evolution of the system. The sites al- 
ready visited by the particle arc connected with the cur- 
rent particle location by an arrow path and the cluster 
of these sites has an acyclic structure. The cluster of 
acyclic arrows is embedded into the media of randomly 
distributed arrows. When the particle enters the cluster 
of visited sites it behaves like in the critical state. Each 
time the linear size of the visited cluster increases by A, 
the particle visits all sites of the cluster four times again. 
The time of increasing A is of order of the size of the 
cluster, i.e. 



dT 
dR 



(18) 



Thus, instead of the simple diffusion law (|T^ in the crit- 
ical state, for the transient states we obtain 



(19) 



which has already been determined in |^ . Note that the 
power law (noh is valid only at the time scale much greater 



than the time spent inside one cluster of the visited sites. 
Inside the cluster, the motion of the particles is similar 
to that in the critical state with the diffusion law (|T^). 
By averaging over a large number of returns to the origin 
one obtains vt — \- 

Now we can estimate the average time required to 
reach the critical state starting from an arbitrary ran- 
dom configuration of arrows. In order to get a spanning 
tree on the lattice, the particle must visit all sites at least 
once. Using (|l^) we can obtain for the lattice of the size 
Lx L 



(20) 



The same time is required for the particle walking on the 
closed graph to settle into the Eulerian trail. 

We also measured numerically the motion of the parti- 
cle in the system. Starting from the transient state, the 
mean square displacement of the particle is described by 
the power law with the critical exponent ut = 0.33 as is 
shown in Fig. pi,. The subsequent evolution of the system 
by the repetitive addition of particles changes this power 
law. For the system in the SOC state, we obtained the 
value V = 0.5 (Fig. ^). These simulations illustrate very 
well the exact results obtained above. 

In summary, we considered the dynamics of the Eule- 
rian Walkers Model. The structure of avalanches in the 
SOC state was studied in detail. We obtained the criti- 
cal exponent for the distribution of durations of the first 
avalanche. Considering the evolution of the system as a 
sequence of avalanches, we found the simple diffusion law 
for the mean square displacement of the particle in the 
SOC state. The crossover from the transient state into 
the SOC state was described as well. The obtained exact 
results were confirmed by numerical simulations. 
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FIG. 1. (a) - closing the loop at ii. (b) - the last step before openning the loop, (c), (d), (e) - the evolution on the closed 
graph settled into the Eulerian trail. The loops in (a) and (b) exactly coincide with those in (d) and (e). 

FIG. 2. The distribution of duration of the first avalanche in the SOC state is shown on the double logarithmic plot. The 
distribution splits into two parts as described in the text. The slope of both the parts is the same with the critical exponent 
T = 2.0. 

FIG. 3. A subsequent evolution of a cluster of visited sites in the SOC state. A schematic picture of the cluster after the first 
(a) and second (b) stages of evolution. The areas with different numbers of visits are shown by different colors. The directions 
of arrows correspond to their final positions. 

FIG. 4. The dependence of the mean square displacement of the particle on time in the transient (a) and SOC (b) states. 
The obtained values of the critical exponents are vt = 0.33 and v = 0.5, respectively. 
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Fig. 4a 
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Fig. 4b 



